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Abstract 

When the bilocal heavy-quark effective theory (HQET) operator for the S-meson 
distribution amphtude has a light-like distance t between the quark and antiquark 
fields, the scale ~ 1/t separates the UV and IR regions, which induce the cusp singu- 
larity in radiative corrections and the mixing of multiparticle states in nonperturbative 
corrections, respectively. We treat these notorious UV and IR behaviors simultaneously 
using the operator product expansion, with the local operators of dimension d < 5 and 
radiative corrections at order ag for the corresponding Wilson coefficients. The result is 
derived in the coordinate space, which manifests the Wilson coefficients with Sudakov- 
type double logarithms and the higher-dimensional operators with additional gluons. 
This result yields the i?-meson distribution amplitude for t less than ~ 1 GeV~^, in 
terms of A = — nify and the two additional HQET parameters as matrix elements 
of dimension-5 operators. The impact of these novel HQET parameters on the integral 
relevant to exclusive B decays, Ab, is also discussed. 



For the exclusive S- meson decays, such as B ^ tttt, p7, . . ., systematic methods have been 
developed using QCD factorization based on the heavy-quark limit [1-3]. Among the building 
blocks in the corresponding factorization formula of the decay amplitude, essential roles are 
played by the light-cone distribution amplitudes (LCDAs) for the participating mesons, 
which include nonperturbative long-distance contributions. In particular, in addition to the 
LCDAs for the light mesons tt, p, etc., produced in the final state, those for the B meson [4] 
also participate in processes where large momentum is transferred to the soft spectator quark 
via hard gluon exchange [1-3, 5-8]. The leading quark-antiquark component of the i?-meson 
LCDA is defined as the vacuum-to-meson matrix element [9]*: 

Here, the bilocal operator is built of the 6-quark and light-antiquark fields, /i„(0) and q(tn), 
linked by the Wilson line at a light-like separation tn, with as the light-like vector {rP = 0, 
n ■ V = 1), and representing the 4- velocity of the B meson; a difference between (1) and 
the familiar pion-LCDA [11] is that /i„(0) is an effective field in the heavy-quark effective 
theory (HQET) [12]. fi denotes the scale where the operator is renormalized, and F{fi) 
is the decay constant in HQET, F{fi) = —i{0\qi^'~f5hj,\B{v)). The RHS in (1) defines the 
momentum representation, with ujv~^ denoting the LC component of the momentum of the 
light antiquark. 

The "IR structure" of the S-meson LCDA was studied [13] using constraints from the 
equations of motion (EOM), qp = v ■ Dh^ = 0, and heavy-quark symmetry, '^hy = hy. 
These constraints are solved to give (1) as 4>+{^) = (f)^^^\uj) + (f)^^\uj), where the first term 
is expressed by the analytic formula (j)^^^\uj) = uj6{2A — u})/{2K^) in terms of a HQET 
parameter [12], A = — m;,, which represents the mass difference between the 5-meson 
and 6-quark. The second term, is obtained as a certain integral of the matrix element 
of the three-body LC operator, {q{tn)Gap{un)a\r^hy{Q)) , where^ the nonperturbative gluons 
participate as the field strength tensor Gq,/3; see (9) below for the integro-differential equation, 
which is derived from the relevant EOM constraints and gives (t)+{oo) = (j)^^^\uj) + ^^■'(co') 
as its solution (see also [14]). In addition, the "UV structure" of the i?-meson DA was 
studied [15] by calculating the 1-loop renormalization of the bilocal operator in (1). The 
"vertex-type" correction around a "cusp" between the two Wilson lines, the light-like Wilson 
line of (1) and the time-like Wilson line from hy{0) = P exp[ig J_^dXv ■ A{Xv)]h^{—oov), 
produces the "radiation tail", given by (f)^{uj,fi) ~ — Cira^ ln(ci;//i)/(7rci;) for u ^ fi, where 
Cp = {N^ — l)/{2Nc). This implies that the moments duu^cp^ioj, jj), which would 
correspond to matrix elements of the local operators q{Q){n ■ Dy'i^'y^hy{0), are divergent, 
reflecting "cusp singularity" [15,16]. 

However, these results do not represent the whole story. As we shall demonstrate ex- 
plicitly, a full description of (1) actually involves a complicated mixture of the IR and UV 
structures; e.g., the above-mentioned functional forms (j)^^^\uj), (f)^^\uj) from the EOM 



*For three-quark LCDAs for the Af, baryon, see [10]. 

tPor simphcity the Wilson hnes connecting the constituent fields are suppressed, and (•■•) = (0| ' l^i"^))- 
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constraints are subject to additional effects from radiative corrections when combined with 
the UV structure and are profoundly modified, in particular, for large values of u (see the 
discussion below (9)). To investigate the behavior of the 5-meson DA incorporating both 
IR and UV structures, we first calculate the radiative corrections, taking into account hard 
and soft /coUinear loops. The one-particle-irreducible 1-loop diagrams (ILDs) for the 2-point 
function {q{tn)'i/i'^^hv(fi)) of (1) yield"'" 



mtn)i-i^hM) - t ( — + 2L - 1 - e ) mtn)v ■ Di^i^hM) 



(2) 



'm D = A — 2e dimensions and Feynman gauge, where L = In [i{t — iQ)ne""^] with the MS 
scale /i and the Euler constant 7e, and the "— iO" prescription comes from the position of the 
pole in the relevant propagators in the coordinate-space (see (3) below). The "vertex-type" 
correction that connects the light-like Wilson line and the q{tn) field in (1) is associated 
with only the massless degrees of freedom; thus, the correction yields the term with the 
"canceling" U"V and IR poles, l/suv — l/^^/R, from the scaleless loop-integral, and with 
the "plus" -distribution (^/(l — 0)+ the splitting function. This term is identical to the 
corresponding correction for the case of the pion LCDA, where the contributions from all 
1-loop diagrams have the same l/suv ~ ^/^ir structure with the Brodsky-Lepage kernel [11] 
as the total splitting function. However, the other terms in (2) have "non-canceling" UV and 
IR poles: another vertex-type correction, connecting the light-like Wilson line and the hy{0) 
field in (1), gives [15] the terms proportional to 5(1 —C,), which contain the double as well as 
single UV pole, corresponding to the cusp singularity mentioned above. The "ladder-type" 
correction, connecting the two quark fields in (1), gives all the remaining terms in (2), which 
contain the IR poles and are associated with not only the bilocal operator in (1), but also 
the higher dimensional operators (see the discussion below (3)); note that the ellipses in (2) 
are expressed by the operators involving two or more additional covariant derivatives. 

The renormalized LCDA is obtained by subtracting the UV poles from (2) with the 
trivial quark self-energy corrections complemented. Here, the term with the plus-distribution 
{^/{l — 0)+ is analytic (Taylor expandable) at t = 0, similar to the pion LCDA, but the 
other terms are not analytic due to the presence of logarithms L, [15, 17]. In particular, 
the nontrivial dependence of the latter terms on tfi through L implies that the scale ~ 1/t 
separates the UV and IR regions. Thus, we have to use the operator product expansion 
(OPE) to treat the different UV and IR behaviors simultaneously: the coefficient functions 
absorb all the singular logarithms, while, for the local operators to absorb the IR poles, 
we have to take into account many higher dimensional operators. Such OPE with local 
operators is useful when the separation t is less than the typical distance scale of quantum 
fluctuation, i.e., when t < The OPE result can be evolved from a low /i to a higher 

scale using the Brodsky-Lepage-type kernel and the Sudakov-type operator with the cusp 
anomalous dimension [15, 17] associated with the single-UV-pole terms in (2). 
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An OPE for the S-meson LCDA (1) was discussed in [18], taking into account the local 
operators of dimension d < A and the NLO (0(as)) corrections to the corresponding Wilson 
coefficients in a "cutoff scheme", where an additional momentum cutoff Aj/y ^qcd) was 
introduced, and the OPE, in powers of l/Af/y, was derived for the regularized moments, 
Jq^^ du!U!^ (I)^{uj , fi) , in particular, for the first two moments with j = 0, 1. In this Letter, 
we derive the OPE for (1), taking into account the local operators of dimension d < 5 and 
calculating the corresponding Wilson coefficients at NLO accuracy. Following the discussion 
above, we carry out the calculation for t < in the coordinate space and in the MS 
scheme, so that there is no need to introduce any additional cutoff. 

The most complicated part of this calculation is the reorganization of contributions from 
(many) Feynman diagrams in terms of the matrix element of gauge-invariant operators, 
including higher dimensional operators. In particular, to derive the Wilson coefficients as- 
sociated with the dimension-5 operators such as qGapi^jbh^, we have to compute the 1-loop 
diagrams for the 3-point function, as well as those for the 2-point function as in (2), where 
the former diagrams are obtained by attaching the external gluon line to the latter dia- 
grams in all possible ways. To perform the calculation in a systematic and economic way, 
we employ the background field method [19]. We decompose the quark and gluon fields into 
the "quantum" and "classical" parts, where the latter part represents the nonperturbative 
long-distance degrees of freedom inside the B meson as a background field and satisfies the 
classical EOM. The quark and gluon propagators for the quantum part contain their cou- 
pling with an arbitrary number of background fields, and each building block of the Feynman 
diagrams obeys the exact transformation property under the gauge transformation for the 
background fields. We use the Fock-Schwinger gauge, x'^A^i^\x) = 0, for the background 
gluon field This gauge condition is solved to give aI^\x) = duux^G^^^^{ux) [19], 

which allows us to reexpress each Feynman diagram in terms of the matrix element of the 
operators associated with the field strength tensor. Using this framework, the tree-level 
matching to derive our OPE can be performed replacing each constituent field in (1) with 
the corresponding background field. The classical bilocal operator can be Taylor expanded, 
and we obtain the OPE at the tree level with the local operators of dimension-3, -4, and -5 
(see the 0{a^) terms in (4) below). 

For the 1-loop matching, we calculate the 1-loop corrections to the 2- and 3-point func- 
tions with the insertion of the bilocal operator in (1), taking into account the mixing of 
operators of dimension d < 5. Apparently, the mixing through the UV region of the loop 
momenta arises only in the 2-point function, and the result can be immediately obtained 
from (2); however, additional mixing can arise in both 2- and 3-point functions accompa- 
nying the IR poles. We perform the loop integration in the coordinate space using the 
Schwinger ( "heat kernel" ) representation of the Feynman amplitudes under the background 
fields [19]. For the calculation of the 3-point function, the Fock-Schwinger gauge ensures 
that the Wilson line in (1), as well as the heavy-quark field, does not couple directly to 
the background gluons, while a massless quark or gluon line couples to them. For exam- 
ple, the "ladder-type" correction diagrams for the 2- and 3-point functions are obtained by 
connecting the two quark fields in (1) using the gluon propagator in the background gluon 
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fields [19], and, after some manipulations using the heat-kernel expansion, we encounter the 
(coordinate-space) loop integrations of the type 



dr- 
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in D (= 4 — 2e) dimensions, with n,m = 1, 2, . . ., and Sn = X]fc=i(-'-/^)- Here, the integral 
over the proper time r, associated with the propagation of the heavy quark in the corre- 
sponding diagrams, is UV finite but IR divergent as -D — > 4, yielding the RHS. The 2-point 



function receives the contribution of the dimension-3 operator. 



,hy, accompanying (3) 



with n = m = 1 as its coefficient. The higher dimensional operators with additional co- 
variant derivatives and/or the gluon field strength tensor contribute to 2- as well as 3-point 
functions, accompanying (3) with n + m > 2. Those latter contributions manifest the mix- 
ing of operators of dimension-4 and higher, accompanying the IR poles. Note that this type 
of mixing mechanism for "higher twist" operators through loop effects is forbidden if the 
corresponding loop subdiagram is composed of massless degrees of freedom only. 

After working out the loop integral of all 1-loop diagrams for the 2- and 3-point functions, 
we reorganize the result in terms of a complete set of gauge-invariant operators of dimension 
d < 5, using some extension of the technique of [9, 13] based on the EOM and heavy-quark 
symmetry. We subtract the UV poles to renormalize the bilocal operator of (1) and also 
renormalize the local operators to absorb the IR poles. (The details of the above-mentioned 
procedures to treat the 2- and 3-point functions at 1-loop will be presented elsewhere [20].) 
Combining the result with the tree-level result obtained above, we obtain the OPE for 
the bilocal operator in (1), g(tn)Pe*3/o '^^"•^(^")^75/i„(0) = ^. Ci(t, /i)C>i(/i), to the desired 
accuracy as 
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Here and below, Cq = N^, fi is the MS scale, and = Q;s(/i). A basis of local operators of 
dimension-(i, Cf (/c = 1, 2, . . .), is defined as 
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(4) 



q{in ■ D)i/L-f5hy , q{w ■ D)ii-i^h^ 



(5) 
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{of} = { 



q{in ■ D^'^'j^h^ , q{iv ■ D){in ■ D)'ifi'~^^h^ , q{iv ■ D] 



qigGaf3V°'n^'^'y5hv 



qigGapl^n^ 



'^'j^hv , qigGapTv'^T/i'-ir^h^ , qgGaf3(T°''^'^l5hv} ,(6) 



with another hght-hke vector, = 0, as = (n^ + n^)/2. The double logarithm in the 
coefficient functions originates from the cusp singularity (see (2)). The 1-loop corrections 
for the 3-point function induce only O4 5 g 7 associated with the field-strength tensor, while 
those for the 2-point function induce all ten operators of (5), (6) through the use of the 
EOM. 

Taking the matrix element (■ ■ ■ ) = {0\ ■ ■ ■ \B{v)) of (4), we can derive the OPE form of 
the S-meson LCDA (1). The matrix elements of the local operators in (5), (6) are known to 
be related to a few nonperturbative parameters in the HQET, using the EOM and heavy- 
quark symmetry as demonstrated in [9, 13]: {0[^^) = 4zF(/i)A/3, {O^^) = iF{jj)K, where F 
and A were introduced below (1), and all seven matrix elements (0^^^) for the dimension-5 
operators (6) can be expressed by F, A and two additional HQET parameters A^; and A//, 
which are associated with the chromoelectric and chromomagnetic fields inside the B meson 
as {qgE ■ (xj^hy) = F(/i)A|;(/i) and {qgH ■ a^y^hy) = iF{fi)X'j^{fi), respectively, in the rest 
frame where f = (1, 0). As a result, we obtain the OPE form for the LCDA (1), 
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which takes into account the Wilson coefficients to 0{as) and a complete set of the local op- 
erators of dimension c? < 5; the cusp singularity in the former leads to the double logarithms 
L^, while the constraints on matrix elements of the latter from the EOM and heavy-quark 
symmetry allow us to represent the result completely in terms of only three HQET param- 
eters, A, A^; and \h- Thus, (7) "merges" the UV [15] and IR structures [13] peculiar to 
the 5-meson LCDA, so that it embodies novel behaviors that are completely different from 
those of the pion LCDA: yU and t are strongly correlated due to the logarithmic contributions, 
L = In [z(t — iO)fie""^], from radiative corrections, so that the DA is not Taylor expandable 
about t = 0. The DA receives the contributions from (many) higher dimensional operators, 
in particular, from those associated with the long-distance gluon fields inside the 5-meson. 

Fourier transforming (7), we obtain the momentum representation, (j)^{uj,fi) in (1), and 
can also evaluate the regularized-moments, Mj = J^^^ duu^ (j)^{uj , fi) , with the cutoff Auv 
(^ Aqcd); Mj corresponds to the "regularized coefficients" of the P term in the Taylor 
expansion of the DA (1) and thus represents the P Mj = O (Mj/A^y) effects on 0+(t, /i) for 
t < 1/A;7y. Here, without going into the detail [20], we note that, for the ffist two moments 
Mq and Ml, the contributions from the ffist line in (7), associated with the matrix elements 
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of the dimension-3 and -4 operators (5), coincide with the result obtained in [18], while the 
second and third lines"'" in (7), which are generated from the dimension-5 operators (6), yield 
the new power-correction terms down by Aqcd/Ac/v, as controlled by the corresponding 
Wilson coefficients. A similar pattern is observed for the asymptotic behavior of 
for uj S> Aqcd- In view of our coordinate-space OPE approach, the UV divergence in the 
moments [9,15,18], Mj — > cxd as Kuv oo, reflects the non-analyticity of (4) and (7) at 
t = due to the presence of logarithms L and in the Wilson coefficients (see also the 
discussion in [17]). 

Our OPE results (4) and (7) give a model-independent description of the S-meson LCDA 
when t < l//i (< I/Aqcd), taking into account the UV and IR structures simultaneously. It 
is instructive to draw a comparison with the previous results mentioned below (1), concerning 
UV [15] or IR structure [13]. In [15], the renormalization group (RG) evolution of (1) was 
derived. The corresponding evolution kernel is determined by the (single) UV poles in (2) as 
the sum of the cusp anomalous dimension and the plus-distribution terms. One can prove 
that our LCDA (7) satisfies the RG equation of [15], taking into account dA/d^ = and [21] 

^d^^\\l{^l)) 47rV|Ci.-|CG ICf + ICg)\\1{^) 

Thus, our results, (4) and (7) via the 1-loop matching, are completely consistent with the 
UV structure obtained at the 1-loop level in [15]. As a result, the scale dependence of the B- 
meson LCDA (1) with t < 1/Aqcd is represented by the two-step evolution: it is governed by 
the solution of the evolution equation in terms of the Sudakov-type and Brodsky-Lepage-type 
scale dependence from the high scale /Xj ~ ^/^^MA^jcD^ associated with the QCD factorization 
formula, to the scale /i < 1/t, while that for the lower scale is governed by the anomalous 
dimensions of local operators like (8). Here, in principle, we can take advantage of the RG 
improvement as usual. To achieve the control at the next-to-leading logarithmic accuracy, 
we take into account the 2-loop as well as 1-loop cusp anomalous dimension [16] for the 
Sudakov-type evolution from Hi to /i (see, e.g., [6]); similarly, the evolution for the lower 
scale requires the 2-loop anomalous dimensions of \\ and A|^, but they are unknown at 
present. 

On the other hand, the connection of the present result (7) with the IR structure of 
[13] is not straightforward: from the EOM and heavy-quark symmetry constraints, a set 
of relations between the two- and three-particle LCDAs was obtained in [13]. Using these 
relations, we can derive the following integro-differential equation for the LCDA (1), 

^<^4f)_ ^ ^2iAt + 1) 0+(t) - - / dt'4>+{t') + 2t^ [ dul (u + l) ^A{t, u) 
dt t Jq jo I 



+ u^v{t, u) + XA{t, ^) - ^ ^ dt't'^ [^A{t', u) - ^vit', m)] I = , (9) 



■'"They also allow us to derive M2, up to the corrections that are suppressed for Auv ^ Aqcd [20]; from di- 



mension counting [18], the local operators of dimension d (> 3) contribute to Alj as ~ A^y(AQCD/A, 
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where and XA{t,u) are the three-particle LCDAs introduced in [13] (see 

also [14]): when t < 1/Aqcd, we substitute (7) for 0+, while the three-particle LCDAs are 
given as = l-^l;) "^v = I-^h; ^^'^ -^a = in terms of matrix elements of local operators 
of dimension d < 5 [13], omitting the contributions from higher dimensional operators as 
well as from radiative corrections. Then the contribution to the LHS of (9) vanishes 

up to the corrections of O ((tAqcD)'^) which are beyond the present accuracy, but the 0{as) 
contribution proves to yield a nonzero result as {asCF/4:7i) [— 8L + it A (32L/3 -|- 3/2) -|- ■ ■ ■ ] . 
Thus, (9) receives corrections at order a^; namely, the relations of [13] are satisfied by the 
nonperturbative matrix elements of local operators in the OPE (4), but the 0{as) loop 
effects in the corresponding Wilson coefficients prevent those relations from being exact for 
the DA (7) (see also the discussion in [17,22]). Such a violation of relations of the type (9) 
at order as in perturbation theory is peculiar to the heavy-meson LCDAs in the HQET and 
does not arise for the case of the (higher twist) LCDAs for the light mesons, vr, p, etc. [23]. 

Our OPE form (7) allows us to parameterize all nonperturbative contributions in the 
5-meson LCDA (1) for t < 1/Aqcd by a few HQET parameters. Here, we evaluate (7) at 
the scale p = 1 GeV, using the information available for those HQET parameters. First of 
all, we note that A = — rrif, in (7) is defined by the 6-quark pole mass m^. Following 
[18], we eliminate A in favor of a short-distance parameter, Ada, free from IR renormalon 
ambiguities [24] and written as A = AoA^fJ') [1 + {7/ 1611)0 pas] — {9/ 811)^10 pas, to one-loop 
accuracy; AnAifJ') can be related to another short- distance mass parameter whose value 
is extracted from analysis of the spectra in inclusive decays B Xg'j and B — »• X^l v, 
leading to A£)^(/i = 1 GeV) ~ 0.52 GeV [18]. For the novel parameters associated with the 
dimension-5 operators, we use the central values of 

A|(/i) = 0.11 ±0.06 GeV^ , A|(/i) = 0.18 ± 0.07 GeV^ , (10) 

at /i = 1 GeV, which were obtained by QCD sum rules [9]; no other estimate exists for 
A^; or \h- We calculate (7) for imaginary LC separation, performing the Wick rotation 
t -> -IT [9,17]. 

The results for 0+(— ir, /i = 1 GeV) using (7) are shown as a function of r in Fig. 1: 
the RHS shows the behaviors for small and moderate r regions, and the LHS displays also 
the region with larger r. The wide-solid curve shows the whole contributions of (7), while 
the narrow-solid curve shows the result for — > 0; the NLO perturbative corrections are 
at the 10-30% level for moderate r, while they are very large for r — > because of singular 
logarithms and L. The dashed and dot-dashed curves show the contributions of the first 
two terms and the first line in (7), respectively, associated with the operators of dimension 
d = ?> and < 4, while the dotted curve gives the results of (7) when Xe = \h = ^- For 
moderate r, the contributions from the dimension-4 operators suppress the DA by 30-40%, 
but the dimension-5 operators, in contrast, lead to enhancement by 10-20% with significant 
effects from A^ and \h- 

In the large r region, the dimension-4 and -5 operators change their roles such that they 
enhance and suppress the DA, respectively, because of the growth of the double logarithm 
in the NLO Wilson coefficients; this fact would eventually call for treating the higher order 
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Figure 1: The 5-meson LCDA (j)^{—iT, fi) at /i = 1 GeV using the OPE-based form (7): 
the wide-sohd curve represents the whole contributions, and the dashed and dot-dashed 
curves show the contributions associated with the operators of dimension d = 3 and d < 4, 
respectively. The narrow-solid and dotted curves show the cases when — > and Xe,h ~^ 0, 
respectively. The two-dot-dashed curve plots the result using the Lee-Neubert ansatz (11). 



perturbative corrections beyond NLO, in particular, resummation of the large logarithms in 
the Wilson coefficients. We also note that the contributions from the dimension-d operators 
grow as ~ t'^~^ with increasing r, and, beyond a certain value of r, the contributions from the 
dimension-4 or -5 operators become larger than those from the lower-dimensional operators; 
this indicates that the hierarchy of contributions is lost, and the OPE breaks down. These 
considerations with quantitative results in Fig. 1 show that our S-meson LCDA (7) indeed 
works up to moderate LC distances r of order 1 GeV^^ ~ l//i; in this region, the hierarchy 
among the dashed, dot-dashed, and wide-solid curves demonstrates convergence of the OPE 
(4) and the corresponding accuracy of our LCDA (7). 

We find that the behavior of the wide-solid curve for small and moderate r in Fig. 1 is 
similar to that of the previous results [17, 18] for the 5-meson LCDA, obtained by taking into 
account perturbative as well as nonperturbative QCD corrections in a systematic framework. 
For example, the two-dot-dashed curve in Fig. 1 shows the behavior of the two-component 
ansatz by Lee and Neubert [18], which is given in momentum space as 



(l)f{u, l^) = N^ e^'^Z"" + 9{iu - tut) 
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where the second term reproduces the correct asymptotic behavior of the DA (1) for u ^ 
Aqcd, and the first term represents the nonperturbative component modeled by an expo- 
nential form [9]; uUf is chosen such that (11) is continuous. The other parameters, and 
ujq, are fixed by matching the first two (j = 0,1) cut-moments J^^^ dujuu^ (f)^ {lu , fi) with 
the OPE for the corresponding cut-moments Mq^i derived in [18], where the operators of 
dimension d < 4 and the corresponding Wilson coefficients at NLO are taken into account 
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(see the discussion below (7)). The (central) values of these parameters are ut = 2.33 GeV, 
= 0.963, and = 0.438 GeV at /i = 1 GeV [18]. For r < 1 GeV^^ in Fig. 1, the Lee- 
Neubert ansatz (11) shows behavior similar to (7) with Xe = Xr = ^ substituted. Indeed, 
when Fourier transformed to the coordinate space, (11) produces, with good accuracy, the 
terms in the first line of (7), which are associated with the operators of dimension d < 4; 
moreover, the first term of (11) also produces particular contributions associated with the 
operators of dimension-5 and higher (see (12) below), and the sizes of those contributions 
are actually rather close to those of the terms proportional to in (7) for r < 1 GeV~^. 

The wide-solid curve in Fig. 1 represents the model-independent behavior of the B- 
meson LCDA (1) based on the most accurate OPE (4), (7) at present. However, in the 
long-distance region, r ^ 1 GeV~^, the contributions associated with the operators of any 
higher dimension become important, and the OPE diverges; thus, one has to rely on a 
certain model for the large r behavior and connect the model-independent descriptions at 
small and moderate r to that model in a reasonable manner. The results in Fig. 1 suggest 
the possibility of connecting the behavior for t < Tc {tc ^ 1 GeV~^) given by our OPE form 
(7) to that for r > Tc, given by the coordinate-space representation of the first term of (11), 

due-''- N ^ e-"/"° = 2 • (12) 







{tujo + 1) 



Here, N and can be determined such that both the resulting total DA (j)^{—iT, jj) and its 
derivative 90+(— zr, /i)/9r are continuous at r = r^- Namely, we perform the matching of 
0+(— zTcyu) = (iu;e~'^'^'=0+(u;, /i), as well as of (90+ (—zTc, /i)/9rc = — duje~'''^''uj(f)j^{uj, ii), 
between (7) and (12), and this is formally analogous to the matching used for (11). In the 
LHS of Table 1, associated with the central values of (10), we show the values of A^ and ujq 
obtained by solving our matching relations for /i = 1 GeV. (The RHS of Table 1 shows the 
results that would be obtained by solving the similar matching relations with A^; = Xh = 0, 
and we find that, indeed, Tc ^ 0.7 GeV~^ gives the behavior close to (11).) In Fig. 2, the wide- 
solid and two-dot-dashed curves are same as those in Fig. 1, and the dotted, solid-gray, and 
dashed curves show the behavior of (12) for r > with t^. = 0.6, 1.0, and 1.4 GeV"^, respec- 
tively, using the corresponding values of A^ and Uq in the LHS of Table 1; these three curves 
behave as ~ N/ {uJqt'^) at large r, with larger N/uJq than those of (11) and the RHS of Table 1. 
Indeed, we can show that N/ul = {9/4Al^) {l + t.Ada [)^1/^Ia + ^^h/C^^Ia) - l] } + ■ ■ ■ , 
using our matching relations, and thus the contributions of A^; and Xh enhance N/uJq. 
Using these results, we calculate the first inverse moment of the LCDA, 

^B^if^) = / '^^ ' — = / dT(p^{—iT, fi) + / dT(f)+{—iT, n) , (13) 

Jo ^ Jo Jtc 

which is of particular interest for the QCD description of exclusive S-meson decays. We 
substitute (7) and (12) into the first and the second terms in the RHS, respectively, and 
the results are shown in Table 1 for each value of Tc- The "stable" behavior observed for 
0.6 GeV"^ < < 1 GeV"^ in the LHS of Table 1 and in Fig. 2 suggests that Xb^{^i = 
1 GeV) ~ 2.7 GeV~\ i.e., Aij(/i = 1 GeV) - 0.37 GeV. This value of Xb is somewhat 
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0.4 
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3.11 (0.23 + 2.88) 
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0.832 


0.301 


2.69 (0.23 + 2.46) 
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U.oUD 
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u.oyy 


n QO/i 
U.oy4 


z.iy (^u.oo + i.o4j 


0.8 


0.852 


0.308 


2.69 (0.47 + 2.22) 


0.966 


0.461 


1.99 (0.46 + 1.53) 


1.0 


0.858 


0.313 


2.66 (0.58 + 2.08) 


1.11 


0.572 


1.79 (0.56 + 1.23) 


1.2 


0.910 


0.349 


2.51 (0.67+ 1.84) 


1.55 


0.839 


1.56 (0.64 + 0.92) 


1.4 


1.09 


0.456 


2.22 (0.76 + 1.46) 


4.43 


1.95 


1.32 (0.71 + 0.61) 


1.6 


1.81 


0.777 


1.87 (0.83 + 1.04) 


9.82 


-4.55 


1.11 (0.77 + 0.34) 



Table 1: Parameters of the model function (12) for different values of Tc, determined by 
continuity at r = with the OPE-based LCDA (7) for yU = 1 GeV, and the results of the 
inverse moment A^^(/i) at /i = 1 GeV, with the first and second numbers in the parentheses 
denoting the contributions from the first and the second terms in the RHS of (13). 

0.8 



> 



=^0.4 

M 
■H 

0.2 

+ 



1 2 3 4 5 

z [GeV"^] 




Figure 2: The S-meson LCDA (j)^{—iT, fi) at /i = 1 GeV. The solid-black and two-dot- 
dashed curves show the results using the OPE-based form (7) and the Lee-Neubert ansatz 
(11), respectively. The dotted, solid-gray and dashed curves plot the model function (12) 
for T > Tc with Tc = 0.6, 1.0 and 1.4 GeV~^, respectively, determined by continuity with (7). 

smaller than the previous estimates that include nonperturbative and/or perturbative QCD 
corrections [17,18,25] (e.g., (11) gives XsifJ' = 1 GeV) ^ 0.48 GeV), although consistent 
with them within their theoretical errors. A value of A^ that is as small as our value was 
adopted in [1,7]. Note that in the RHS of Table 1 with Xe,h = 0, the stable behavior is not 
seen as clearly as in the LHS, and A^ assumes larger values than in the latter. 

The above results demonstrate that the novel HQET parameters, A^; and Xh, associated 
with the dimension-5 quark- antiquark-gluon operators, could lead to smaller value of A^. In 
particular, this effect can be more significant: using the values of A^ and Xh corresponding 
to the upper bound of (10), we find that the solid curve in Fig. 2 becomes further enhanced 



10 



in the moderate r region, so that (13) gives Xsif^ = 1 GeV) ~ 0.2 GeV or smaller. This 
finding calls for more precise nonperturbative estimates of A^; and A^- We also note that in 
the RHS of (13) evaluated in Table 1, the second term is much larger than the first term. 
This suggests that A^ is rather sensitive to the functional form that models the LCDA in the 
long-distance region; for example, a functional form motivated by which is mentioned 

below (1), provides an interesting possible alternative to (12). Moreover, the Sudakov-type 
resummation for the Wilson coefficients may give rise to additional important contribution 
in the large r region. Systematic investigations of these points, as well as the effect of the RG 
evolution of the LCDA, are beyond the scope of this Letter and will be presented elsewhere. 

To summarize, we have derived the OPE that embodies both the notorious UV and IR 
behaviors of the 5-meson LCDA, including all contributions from the local operators of 
dimension d < 5 and the corresponding Wilson coefficients at NLO accuracy. This OPE 
allows us to parameterize all nonperturbative contributions in terms of three HQET param- 
eters and provides us with the most accurate description of the 5-meson LCDA for distances 
less than ~ 1/Aqcd- We have also used the model-independent behaviors from our OPE 
to constrain the long-distance behavior of the LCDA and estimate the first inverse moment 
A^^ as the integral of the LCDA over entire distances. The results demonstrated the im- 
pact of the novel HQET parameters, associated with the matrix elements of the dimension-5 
quark-antiquark-gluon operators. 
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